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Abstract 



We present spectra of charged pions and protons in 0-10% central Au+Au collisions 
at sj SjVTV = 200 GeV at mid-rapidity (y = 0) and forward pseudorapidity (ij = 2.2) 
measured with the BRAHMS experiment at RHIC. The spectra are compared to 
spectra from p+p collisions at the same energy scaled by the number of binary 
collisions. The resulting nuclear modification factors for central Au+Au collisions 
at both y = and r] = 2.2 exhibit suppression for charged pions but not for (anti- 
protons at intermediate px- The p/ir~ ratios have been measured up to px ~ 3 
GeV/c at the two rapidities and the results indicate that a significant fraction of 
the charged hadrons produced at intermediate pt range are (anti-)protons at both 
mid-rapidity and rj = 2.2. 
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1 Introduction 



One of the reasons for studying heavy-ion collisions at high energies is to 
search for the predicted Quark-Gluon Plasma (QGP), a deconfmed state of 
quarks and gluons, and to investigate the properties of this state of matter at 
extremely high energy densities. High px hadrons, primarily produced from 
the fragmentation of the hard-scattered partons, are considered a good probe 
of the QGP [1,2,3]. Due to induced gluon radiation, hard-scattered partons 
will suffer a larger energy loss in a hot dense medium of color charges than 
in color neutral matter. This results in fewer charged hadrons produced at 
moderate to high px] the hadrons are said to be suppressed. Indeed, all four 
experiments at RHIC have observed that high px inclusive hadron yields in 
central Au+Au collisions are suppressed as compared to p+p and d+Au in- 
teractions at mid-rapidity [4,5,6,7]. However, it was also discovered that the 
yields of protons and anti-protons at intermediate px (1.5-4.5 GeV/c) are com- 
parable to those of pions and not suppressed at mid-rapidity as compared to 
elementary nucleon-nucleon collisions [5,8]. These experimental results have 
motivated several suggestions on how hadrons are produced at intermediate 
Px [9,10,11,12], such as the possibility that boosted quarks from a collectively 
expanding QGP recombine to form the final-state hadrons [10,11,12]. Among 
the interesting results from the BRAHMS experiment is that at forward pseu- 
dorapidity i] = 2.2 inclusive negatively charged hadrons are suppressed in both 
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central Au+Au and minimum-bias d+Au collisions [4,13]. This raises the pos- 
sibility that initial-state effects such as gluon saturation may also influence 
hadron production at intermediate pt [14,15,16]. 

To explore the effect of the nuclear medium on intermediate px particle pro- 
duction, we present in this paper the invariant px spectra of charged pions 
and protons measured by the BRAHMS experiment at RHIC up to 3 GeV/c 
in central Au+Au collisions at v /s/v r /v = 200 GeV at both mid-rapidity and 
forward pseudorapidity r\ = 2.2. The spectra are then compared to reference 
data from p+p collisions at the same energy scaled by the number of binary 
collisions (N bin ) by using the nuclear modification factor: 

d 2 N AA /dp T dy 

((N^/^Jd^/dprdy' [ } 

where d 2 N AA j 'dpxdy is the differential yield per event in the nucleus-nucleus 
(A+A) collision, and of% e i an d d 2 af^ el / 'dprdy are the total and differential cross 
section for inelastic p+p collisions, respectively. 



2 Experiment and data analysis 

The BRAHMS experiment [17] consists of event characterization detectors 
and two independent magnetic spectrometers, the mid-rapidity spectrometer 
(MRS) and the Forward Spectrometer (FS), both of which can be rotated 
in the horizontal plane around the beam direction. For the present studies 
the MRS was positioned at 90 degrees and the FS at 12 degrees with respect 
to the beam direction. Collision centrality is determined from the charged 
particle multiplicity measured by multiplicity detectors as described in [18]. 
The trajectories of charged particles are reconstructed in the tracking devices 
(time projection chambers and drift chambers). The resulting straight-line 
track segments in two detectors located on either side of a magnet are then 
matched and the particle momentum is determined from the deflection of the 
track in the magnetic field. The intrinsic momentum resolution of the spec- 
trometers at maximum magnetic field setting is Sp/p = 0.0077p for the MRS 
and 5p/p = 0.0008p for the FS [13], where p is written in units of GeV/c. In 
the MRS charged particles are identified using a time-of-flight wall (TOFW), 
whereas in the FS a time-of-flight wall (H2) and a ring imaging Cherenkov 
(RICH) detector are used for particle identification (PID). To identify charged 
pions and protons using time-of-flight detectors, 2a standard deviation PID 
cuts in the derived m 2 and momentum space are imposed for each species. 
With a timing resolution of otofw ~ 80 ps in the Au+Au runs, protons and 
pions can be well separated from kaons up to momenta 3.2 GeV/c and 2.0 
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GeV/c, respectively. For pions above 2 GeV/c, an asymmetric PID cut is ap- 
plied, i.e. the region where the pion and kaon 2a cuts overlap is excluded for 
PID, and the pion yield in the region is obtained by assuming a symmetric 
PID distribution about the mean pion mass squared value. This allows the 
pion pt spectrum to be extended to 3 GeV/c, at which point the kaon con- 
tamination of pions is estimated to be less than 5% and is accounted for in the 
systematic errors. For the FS PID in the present analysis, H2 is used only for 
the low momentum data. With a timing resolution of <7h2 ~ 90 ps, protons and 
pions can be identified up to 7.1 GeV/c and 4.2 GeV/c, respectively with a 2a 
separation. Above 7.1 GeV/c, an asymmetric PID cut is applied and the pro- 
ton yields in the overlap region are estimated by assuming a symmetric PID 
distribution about the mean proton mass squared value. Between 7.9 GeV/c 
and 9 GeV/c, the Cherenkov threshold for protons, the RICH detector is used 
to determine the kaon contamination of the proton spectrum. At 9 GeV/c the 
contamination of protons by kaons is estimated to be less than 6%. Above 9 
GeV/c, protons are identified by using the RICH to veto pions and kaons. To 
identify pions, the RICH is directly used to separated pions from kaons well 
from momentum of 2.5 GeV/c up to 20 GeV/c. 

The invariant differential yields 4 — tt4 (respectively — d ] N , at forward 
rapidity) were constructed for each spectrometer setting. As discussed in [19] 
the differential yields were corrected for geometrical acceptance, tracking and 
PID inefficiencies, in-flight decay of pions, the effect of absorption and multi- 
ple scattering. The pion contamination by hyperon (A) and neutral Kg decays 
were investigated in [20] and found to be less than 5% in the MRS and 7% 
in the FS, respectively. The contribution to proton spectra by the A decays 
was estimated with a GEANT [21] simulation where an exponential distribu- 
tion in pt with inverse slope taken from the PHENIX and STAR measure- 
ments [22,23] for both (anti-)protons and ( ant i-) lamb das was generated for 
several spectrometer settings. By taking the ratio of A (A) to p(p) yields of 
0.89 (0.95) [22] in 0-10% central Au+Au and 0.45 (0.55) [23] in p+p collisions 
at y/sNN = 200 GeV and assuming a constant behavior in the rapidity interval 
\y\ < 2.2 as indicated by HIJING model [24], it is found that the fraction of 
protons originating from A(A) decays is at a maximum value around 35-40% 
in central Au+Au and 27-30% in p+p collisions and decreases with p T . In the 
following correction for feed-down from the (anti-)lambda decays has been 
applied, whereas the contamination of pions due to weak decays has not been 
corrected but is accounted for in the systematic errors. 



3 Particle spectra 

The top row of Figure 1 shows the px spectra of charged pions (left panel) 
and protons (right panel) at mid-rapidity in 0-10% central Au+Au and p+p 
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p T [GeV/c] p T [GeV/c] 

Fig. 1. Top row: px spectra of charged pions (left panel) and protons (right panel) at 
mid-rapidity in 0-10% central Au+Au and p+p collisions at ^snn = 200 GeV. The 
error bars are statistical only. The systematic errors are estimated to be less than 
15% for pions and 18% for (anti-)protons. For the reference spectrum the systematic 
error is estimated to be less than 19%. For clarity, some spectra are scaled vertically 
as quoted. Bottom row: px spectra for tt~ and p at forward rapidity rj = 2.2 in 0-10% 
central Au+Au and p+p collisions at -^/s/vtv = 200 GeV. The systematic errors are 
estimated to be 14% for pions and 17% for anti-protons. 
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collisions at ^snn = 200 GeV. Also shown in the left panel of the figure is 
the measured spectrum of (n + + 7r~)/2 in p+p collisions, where pions can 
only be identified up to 1.5 GeV/c with TOFW for 2003 p+p runs. We thus 
constructed a reference spectrum shown as a solid line by dividing the neutral 
pion spectrum in p+p collisions measured by PHENIX [25] by the spectrum 
from PYTHIA simulation [26] at the same rapidity range and then multiplying 
the results by the (7r + + 7r~)/2 spectrum from PYTHIA. The spectra of (anti- 
protons in p+p collisions are measured by the BRAHMS spectrometer but to 
a smaller p T coverage compared to those in Au+Au collisions due to a worse 
TOF resolution in the p+p runs. The spectra have been corrected for the 
trigger inefficiency [13] and fitted with an exponential function as shown in 
solid lines with the rapidity density and the inverse slope parameter of 0.101 + 
0.004 (0.098 ± 0.004) and 0.304 ± 0.005 (0.285 ± 0.005 GeV) for proton (anti- 
proton), respectively. The error bars are statistical only. The systematic errors 
in the measured spectra, which come from the uncertainties in the momentum 
determination, the time-of-flight measurements and ring radius reconstruction 
procedures, and the uncertainties in the corrections estimations, are estimated 
to be less than 15% for pions and 18% for (anti-)protons. The systematic error 
in the reconstructed reference spectrum for charged pions is estimated to be 
less than 19%. 

The bottom row of Figure 1 shows the pr spectra for tt~ (left panel) and p 
(right panel) at forward rapidity i] = 2.2 in 0-10% central Au+Au and p+p 
collisions at y / Sjv r /v = 200 GeV. Solid lines are curves fit to the tc~ and p 
spectra in p+p collisions. The spectra are constructed in terms of d 2 N/ dp T drj 
because the rapidity coverages of the FS at 12 degrees for pions and protons 
are different, making a comparison of anti-proton to pion yields difficult. In 
addition, since the Jacobian effect is largest at mid-rapidity and gets rather 
small at larger rapidities at intermediate pr range as we focused on in this 
paper, we expect the conclusions drawn from the spectra expressed in terms 
of d 2 N/dpTdrj should be the same as from those of d 2 N/dpTdy. 



4 Nuclear modification factor 

In Figure 2 the nuclear modification factors for n~ and p are deduced for 0-10% 
central Au+Au collisions at ^^snn — 200 GeV at mid-rapidity (left panel) and 
7] = 2.2 (right panel). Error bars represent statistical errors; the systematic 
errors are indicated by horizontal lines. The dotted and dashed lines indicate 
the expectations of participant scaling and binary scaling, respectively. The 
shaded bars represent the systematic errors associated with the determination 
of these quantities. Systematic errors other than the uncertainties in {Ni, in ) 
determinations are estimated to be 20% except for n~ at mid-rapidity, where 
they are around 24%. 
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Fig. 2. Nuclear modification factors for tt~ and p measured for 0-10% central Au+Au 
collisions at ■ v /s/v r /v = 200 GeV at mid-rapidity (left panel) and r\ = 2.2 (right 
panel). Error bars represent statistical errors; the systematic errors are indicated by 
horizontal lines. The dotted and dashed lines indicate the expectations of participant 
scaling and binary scaling, respectively. The shaded bars represent the systematic 
errors associated with the determination of these quantities. Systematic errors other 
than the uncertainties in (iV^) determinations are estimated to be 20% except for 
7r~ at mid-rapidity, where they are around 24%. 



Similar to the unidentified charged hadrons [4] at both mid-rapidity and for- 
ward rapidity, RauAu for charged pions increases monotonically up to ~ 1.5 
GeV/c and levels off at a value below unity above 1.5 GeV/c indicating that 
charged pions yields are suppressed with respect to p+p collisions at interme- 
diate pt- Furthermore, the ir~ yields at forward rapidity show a similar or even 
stronger suppression, indicating that nuclear effects other than parton energy 
loss (jet quenching) might be contributing to the strong suppression. The sup- 
pression at midrapidity around px ~ 2 GeV/c is smaller (about 30%) than the 
suppression that has been reported for neutral pions [27] and which is seen 
at forward rapidity. This difference can - to a large extent - be attributed to 
the construction of the reference spectrum and has been accounted for by the 
systematical error for tt~ at mid-rapidity. Another interesting feature shown 
in the figure is that the anti-proton yields at both mid-rapidity and rj = 2.2 
are not suppressed at px > 1.5 GeV/c. 



7 



I Q. 



0.9 




O y = 


(Au+Au) 


0.8 




□ Tl = 

y- 


2.2 (Au+Au) 
(p+p) 


0.7 




n- 


2.2 (p + p) 



0.6 
0.5 F- 
0.4 
0.3 
0.2 
0.1 

a 



-o- 




o- 



-□- 



1.5 2 

p T [GeV/c] 



2.5 



3.5 



Fig. 3. p/tt~ ratios at both mid-rapidity and rj = 2.2 for 0-10% central Au+Au 
collisions at ^Jsnn = 200 GeV. The error bars show the statistical errors only. The 
systematic errors are estimated to be less than 12% at both y = and rj = 2.2. The 
corresponding ratios in p + p collisions at y / J/v r /v = 200 GeV are sketched as solid 
line and dotted line, respectively. 



5 Particle ratios 



Figure 3 shows p/ir~ ratios at both mid-rapidity and t] = 2.2 for 0-10% central 
Au+Au collisions at ^/snn = 200 GeV. The error bars show the statistical er- 
rors only. In the present ratios, most systematic errors cancel out. Remaining 
systematic errors arising from PID efficiencies, acceptance corrections, correc- 
tions for nuclear interactions with detector etc. are estimated to be less than 
12% at both y = and rj = 2.2. Also shown in the figure are the corresponding 
ratios for p+p collisions at a/sJvtv = 200 GeV. There is a clear increase of the 
P/ti~ ratios at intermediate pt in central Au+Au collisions relative to the level 
seen in p+p collisions (see also [8,28]). This enhancement is most likely due 
to the interplay of several final-state effects and possibly a new hadronization 
mechanism other than parton fragmentation. Calculations based on a parton 
recombination scenario [12] with a collective flow at the partonic level appear 
to be able to qualitatively describe the data at mid-rapidity. 
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6 Summary 



In summary, the BRAHMS measurements demonstrate a significant suppres- 
sion of charged pions at intermediate px at both mid-rapidity and forward 
rapidity for 0-10% central Au+Au collisions at y/SNN = 200 GeV. Such a 
strong suppression is believed to be caused primarily by the parton losing 
energy when traversing the partonic (i.e. characterized by color degrees of 
freedom) medium created in central Au+Au collisions. The strong 7r~ sup- 
pression at forward rapidity suggests that the hot dense partonic medium 
may also exist in the forward rapidity region and that there might be other 
nuclear effects such as gluon saturation contributing to the suppression. How- 
ever, the suppression is not observed for (anti-)protons at intermediate px at 
either mid-rapidity or forward pseudorapidity r\ = 2.2. p/ti~ ratios in cen- 
tral Au+Au show an enhancement of (anti-)proton production relative to the 
p+p collisions at intermediate px- All these observations are consistent with 
a picture where a dense strongly interacting partonic matter with a strong 
collective flow is most likely formed in central Au+Au collisions over a large 
rapidity range which results in the strong suppression of charged pion yields 
and boosts the protons to higher transverse momentum. 

This work was supported by the Division of Nuclear Physics of the Office of 
Science of the U.S. DOE under contract DE-AC02-98-CG10886, the Danish 
Natural Science Research Council, the Research Council of Norway, the Pol- 
ish State Commission for Scientific Research and the Romanian Ministry of 
Research. 
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